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Abstract. We build a sample of O VI absorption sytems in the redshift range 2.0 < z < 2.6 
using high spectral resolution data of ten quasars from the VLT- UVES Large Programme. We 
investigate the existence of a metal-rich O VI population and define observational criteria for this 
class of aborbers under the assumption of photoionization. The low temperatures of nearly half 
of all O VI aborbers, implied by their line widths, are too low for collisional ionization to be a 
dominant process. We estimate the oxygen abundance under the assumption of photoionization; 
a striking result is the bimodal distribution of [O/H] with median values close to 0.01 and 
0.5 solar for the metal-poor and metal-rich populations, respectively. Using the line widths to 
fix the temperature or assuming a constant, low gas density does not drastically change the 
metallicities of the metal-rich population. We present the first estimate of the Ovi column 
density distribution. Assuming a single power-law distribution, /(N) oc N~ Q , yields a ~ 1.7 and 
a normalization of /(N) = 2.3 x 10~ 13 at log N(0 vi) ~ 13.5, both with a ~30% uncertainty. The 
value of a is similar to that found for C IV surveys, whereas the normalization factor is about 
ten times higher. We use /(N) to derive the number density per unit z and cosmic density, 
fib (Ovi), selecting a limited column density range not strongly affected by incompleteness or 
sample variance. Comparing our results with those obtained at z ~ 0.1 for a similar range 
of column densities implies some decline of dn/dz with z. The cosmic Ovi density derived 
from /(N), f2t>(Ovi)~ (3.5±o;|) x 10~ 7 , is 2.3 times higher than the value estimated using the 
observed O vi sample (of which the metal-rich population contributes ~35%) , easing the problem 
of missing metals at high z (~ 1/4 of the produced metals) but not solving it. We find that 
the majority of the metal-rich absorbers are located within ~ 450 km s _1 of strong Ly-c* lines 
and show that, contrary to the metal-poor absorbers, this population cannot be in hydrostatic 
equilibrium. All of the O vi absorber properties imply that there are two distinct populations: 
metal-poor absorbers tracing the intergalactic medium and metal-rich absorbers associated with 
active sites of star formation and most probably linked to galactic winds. 
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1. Introduction 

There are two major open questions that could be solved by the existence of a warm- 
hot and/or highly ionized phase of the intergalactic medium (IGM): the missing baryons 
at low redshift, z ~ 0-0.5, and the missing metals at z ~ 2.5. 

The baryon budget at low z implies that about 45% of the cosmic baryons are still in 
the form of ionized gas in the IGM (Fukugita et al. 1998). The missing baryonic matter 
could reside in a warm-hot IGM (WHIM) as predicted by hierarchical structure formation 
models (see e.g. Cen & Ostriker 1999; Dave et al. 2001) . The cooler phase of the WHIM 
can be probed by O VIAA1031, 1037 absorption but, at the sensitivity of the FUSE and 
HST surveys, the contribution to the cosmic baryon density, £lb , of the detected O VI 
absorbers is only ~5% (Savage et al. 2002; Richter et al. 2004). The hotter phase of 
the WHIM, T > 5 x 10 5 K, can be probed by Ovii-Ovin X-ray absorption. There are 
very few suitably bright targets for X-ray spectroscopy with Chandra and XMM-Newton, 
and the rare confirmed detections, at a significance level higher than 3cr (Nicastro et al. 
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2005a, b), have Ovn column densities about ten times larger than those of Ovi. The 
contribution of this intergalactic hot phase to fib could be up to ten times higher than 
that of the cooler WHIM. 

At z ~ 2.5, at least 90% of the baryons are in the Ly-a forest, but only about 10% of the 
metals produced by star formation activity in Lyman Break Galaxies (LBGs) have been 
detected up to now (Pettini 1999). The mean metal enrichment of the IGM could reach a 
value Z ~ 0.04 Z Q (Pettini 1999) and recent simulations of galactic winds give estimates 
in the range 0.01-0.06 Z Q (Bcrtone et al. 2005). The observed Civ cosmic density equals 
fib (Civ) ~ 7x 10~ 8 (Songaila 2001; Scannapieco et al. 2005) and, assuming an ionization 
correction of about a factor two, leads to a cosmic abundance [C/H]~ —2.9, thus a 
shortfall of metals by a factor of at least ten and maybe up to 10 2 . The missing metals 
could reside in hot gaseous halos around star-forming galaxies (Pettini 1999; Ferrara 
et al. 2005) and the cooler part of these hot bubbles might be traced by O vi absorption. 

A few surveys of Ovi absorbers at z ~ 2.0-2.5 have already been conducted at the 
VLT and Keck telescopes, some for only a limited number of sightlincs (Bergeron et al. 
2002; Carswell et al. 2002; Simcoe et al. 2002,2004). A non- negligible fraction, -1/3, of 
the Ovi absorptions associated with the Ly-a forest have line widths b < 14 km s _1 , 
thus T < 2 x 10 5 K, which favors a radiative ionization process. A hard UV background 
flux, i.e. small discontinuity at 4 Ryd (Haardt & Madau 1996), reproduces well the 
observed ionic ratios for —3.0 < [Z/H] < —0.5. The inferred values of fib(Ovi) of the 
above surveys are « 1.1 x 10 -7 (assuming fi m , fib, h — 0.7, 0.3, 0.04, 70 throughout 
this paper). Applying a conservative ionization correction, Ovi/O=0.15, yields a mean 
oxygen abundance of [0/H]~ —2.7, and thus, as for C iv surveys, leaves open the problem 
of the missing metals. 

However, a higher metallicity has been derived for O v absorbers. The EUV O vA630 
singlet was searched for and detected in a stacked composite absorption spectra from 
HST-FOS data for absorbers at 1.6 < z < 2.9 with a large range of Hi column densities 
(Telfer et al. 2002). Except in the strongest Hi systems, the lack of detection of the 
associated EUV O iv doublet also suggests a hard ionizing background flux, and the 
derived oxygen abundance is [0/H]~ —2.2 to —1.3. 

The paper is organized as follows: our new VLT Ovi sample is presented in § 2. In 
§ 3, we give results on the oxygen abundance derived under various assumptions on the 
ionization process. The O vi column density distribution and the contribution of the O VI 
absorbers to the cosmic baryon density are given in § 4. The origin of these absorbers is 
discussed in § 5. We present our conclusions and prospectives in § 6. 

2. The O vi sample 

The VLT- UVES Large Programme "The Cosmic Evolution of the IGM" provides a 
homogeneous sample of quasar sightlincs, with emphasis given to lower redshift quasars 
(z < 3) to take advantage of the high UV sensitivity of UVES. This allows a study of 
Ovi systems in the range z — 2.0-2.5 where the crowding of the Ly-a forest is not too 
severe. Altogether, the sample comprises 21 bright quasars (most with V < 17), of which 
19 are at 2 < z < 4, observed with dichroics blue and red. The spectral resolution is 
R = 45,000 (line width b = 6.6 km s _1 ) and the exposure time per setting per quasar 
(2 settings per quasar) of 6 to 10 hr yields a signal-to- noise S/N — 30-40 and 100 at 
3300 and 5500 A respectively. The data were reduced using an upgraded version of the 
ESO-UVES data-reduction pipeline (Aracil et al. in preparation). 

We present results derived from the spectral analysis of ten quasars at 2.1 < z em < 
2.8. Our Ovi sample comprises 136 absorbers with column densities in the range 12.7 
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< log N(Ovi) < 14.6. These absorbers span the redshift interval 1.99 < z < 2.57 with 
a mean value of z = 2.28. Due to partial blending of the associated Hi absorptions and 
small velocity differences between the O VI and H I components, we group the O VI and 
Hi absorptions into 51 systems. Absorption systems within 5000 km s _1 of the quasar 
emission redshift are excluded from this sample. 

2.1. The Ovi subsamples 

There are unusual Ovi absorbers with high abundances, —1 < [O/H] < 0, in previous 
Ovi surveys (Bergeron et al. 2002; Carswell et al. 2002). They have high ionic ratios, 
N(Ovi)/N(Hi) > 0.5, and low Hi column densities, log N(Hi) < 13.0. The survey by 
Simcoc et al. (2004) focussed on systems with log N(Hi) > 13.6, which could account 
for the underrepresentation of highly metal-rich Ovi absorbers in their Ovi sample. 
Since these intriguing systems are not present in every sightlinc, a large quasar sample 
is mandatory for a statistically significant number of metal-rich O vi absorbers. 

It is possible to define the class of metal-rich absorbers using observed column density 
ratios derived from photoionization models since the small line widths of a non-negligible 
fraction of the Ovi systems imply "low" gas temperatures (see § 3.1). Adopting a hard 
UV background spectrum together with a 0.1 solar metallicity leads to observational 
identification criteria for the following classes of absorbers: 

- type 1: N(0 vi)/N(Hi) > 0.25: metal-rich absorbers, 

- type 0: N(Ovi)/N(Hi) < 0.25: metal-poor absorbers. 
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Figure 1. Metal column densities of Ovi and Civ versus Hi column density. The dashed and 
dotted line give the locations of systems with N(Ovi)/N(Hi) = 0.25 and N(Crv)/N(Hi) = 
0.015, respectively. 
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Figure 2. Hi and metal absorptions of type 1 Ovi absorbers: a weak Hi system at z = 2.468 
is shown in the left panel, and a strong H I system at z = 2.398 is shown in the right panel. The 
latter is 140 km s _1 away from a type absorber. 
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Figure 3. Hi and metal absorptions of a type Ovi absorber at z = 2.089 (left panel) and a 
type 2 Ovi absorber at z — 2.314 (right panel). 



There are 39 Ovi type 1 components, 12.9 < log N(0 vi) < 14.5, grouped in 14 Ovi-Hl 
systems. 

A similar criterium is derived for the C IV systems and is used to identify C iv-only 
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Figure 4. H I and metal absorptions of C iv-only type 1 absorbers: two low redshift systems at 
z = 1.727 and 1.729 (left panel) and one system at z = 2.415 with Ovi lines fully blended with 
saturated Lyman lines (right panel). 

metal-rich absorbers (Ovi doublet either outside the observing range, z < 2.0, or fully 
blended with saturated Lyman lines): 

- Civ-only type 1: N(Civ)/N(Hi) > 0.015: metal-rich absorbers. 

There are 18 C iv-only type 1 components, 11.8 < log N(Clv) < 13.8, grouped in 8 
C iv-H I systems. 

Finally, a few absorbers with Ovi blended with strong Lyman lines, and thus with 
uncertain values of N(Ovi), are labelled type 2. 

These different classes of absorbers are shown in figure 1. About 70% of the Ovi+ClV- 
only type 1 absorbers have weak associated Hi lines, log N(Hi) < 13.6, and the type 
and type 1 Ovi absorbers span roughly the same N(Ovi) range. This demonstrates the 
importance of searching for Ovi systems whatever the strength of their associated Hi 
aborption. 

Examples of type 1 absorbers are given in figure 2. The proximity in velocity space of 
a strong Ly-a system to a type 1 absorber is investigated in § 5.1. Examples of types & 
2 O VI absorbers and C iv-only type 1 absorbers are given in figures 3 and 4, respectively. 

3. Abundances 

3.1. Ovi line widths 

The histogram of Ovi line widths is shown in figure 5. There are 81, 39 and 16 Ovi 
components for the types 0, 1 and 2 absorbers, respectively. The bulk of the b distributions 
of the types and 1 overlap, but not their high velocity tails. A Kolmogorov-Smirnov 
test shows that these distributions indeed differ at the 98% confidence level. It should be 
stressed that among the broader absorbers, b > 16 km s _1 , most components are blends 
of several Ovi lines within a velocity range of a few tens of km s . Very few individual 
components are unambiguously broad. 
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Close to half (43%) of the O vi absorbers have line widths b < 12 km s _1 . This confirms 
the results previously found with smaller O vi samples. The implied temperatures, T < 
1.4 x 10 5 K, as also found for the non-saturated associated Lyman lines, are too low 
for O VI to be produced by collisional ionization even for abundances close to solar. We 
will thus assume that photoionization is the dominant ionization process, but will also 
consider simple cases where there is additional collisional heating of the gas, possibly 
through shocks. 




Figure 5. Distribution of the Ovi line widths with the totals outlined. 



3.2. Abundances under the assumption of photoionization 

Following previous discussions on constraints of the spectral energy distribution of the 
ionizing background flux (e.g. Bergeron et al. 2002; Carswell et al. 2002; Telfer et al. 
2002), we select a hard UV metagalactic flux (Haardt & Madau 1996) to derive the 
gas ionization level. We used the CLOUDY v94.0 code (Ferland et al. 1998) to estimate 
ionic column density ratios as a function of the ionization parameter, U = n v /nu, and 
assumed solar relative abundances (Anders & Grevesse 1989). For each system, the value 
of U is fixed by the observed ionic ratio N(Ovi)/N(Civ) which is applicable only if 
Ovi and Civ arc in the same phase. This should be the case for most absorbers as 
Si iv is not detected, except in a few systems with large N(Hi) (> 10 15 cm ~ 2 ). The 
observed range of N(Ovi)/N(Civ) implies —1.4 ^ log U < —0.4, thus an ionization 
ratio 0.09 < Ovi/O < 0.21. 

The distributions of the derived oxygen abundances are presented in figure 6 for the 
31, 14 and 6 Ovi systems of type 0, 1 and 2, respectively. Contrary to the distribution 
of 6(0 Vi) shown above, there is very little overlap between the [O/H] distributions of 
the types and 1 populations. This further suggests that they are indeed two distinct 
populations. The type 2 [O/H] distribution spans a small range in between those of the 
other two populations. 

3.3. Abundances under the assumption of other ionization processes 

To confirm the difference in metallicity between the type (population tracing the IGM) 
and type 1 (population tracing highly metal-enriched sites) absorbers, we investigate 
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Figure 6. Distribution of the oxygen abundance in the photoionization case. The [O/H] 
median values for the types 0, 1 and 2 populations are —2.07, —0.33 and —1.56, respectively. 
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Figure 7. Distributions of the oxygen abundance for the type 1 population under various 
ionization conditions: (1) pure photoionization case (top panel), (2) constant gas density (middle 
panel), (3) photoionization plus temparature fixed by &(0 vi) (bottom panel). The [O/H] median 
values for the cases 1, 2 and 3 are —0.33, —0.80 and —0.35, respectively. 



whether other ionization conditions could yield much lower abundances for the type 1 
population. 

First, we consider a gas phase of constant density, p, thus a single value of the ionization 
paramater. We select an overdensity S = (p/p) ~ 10 at z m 2.2 which is within the range 
of values found in the previous O vi surveys. This yields log U — —0.5, thus an ionization 
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ratio Ovi/O = 0.16. In a large fraction of these cases, Ovi and Civ do not trace the 
same phase. 

Secondly, we reconsider photoionization by a hard UV background flux but now with 
the gas temperature derived from the b value of the main O vi component of each system. 
This is to account for possible additional shock heating. The value of U is still derived 
from the observed N(0 Vi) /N(C iv) ionic ratio. There is no solution for absorbers with 
T > 2.0 x 10 5 K, or b > 14 km s _1 , implying that at these higher temperatures Ovi 
and C iv are not co-spatial. However, for the type 1 systems, the main O vi component 
is always narrower than 14 km s" 1 except in the case of one Ovi doublet blended with 
Lyman lines. 

Together with the case of pure photoionization, the distributions of the oxygen abun- 
dances estimated in the above two cases are shown in figure 7. Although the values of 
[O/H] are somewhat lower under the new ionization conditions, they remain far higher 
than those of the type population. This confirms that the types and 1 O vi absorbers 
trace markedly different populations. 

4. Contribution of the O VI absorbers to the cosmic baryon density 

4.1. Column density distribution 

The column density distribution, /(N), of Ovi absorbers per unit redshift path per unit 
column density can be written: 



where n is the number of Ovi absorbers in a column density bin AN centered on N 
for a sample of quasars with total redshift path AX. For our adopted cosmology, the 
redshift path is defined as: 



The Ovi column density distribution is shown in figure 8. It can be seen that the 
present data become incomplete below a column density of ~ 1 x 10 13 cm~ 2 and that 
sample variance may be important at column densities larger than ~ 2 x 10 14 cm~ 2 . 
In between these limits, a power-law fit (/(N) oc N~ Q ) gives a ~ 1.7. To estimate the 
uncertainty in a, we shift the AN bins by 0.1 dex and derive new power-law fits. This 
yields a = 1.71±g;|f and a normalization of /(N) = 2.3 x 10~ 13 at log N(Ovi) = 13.5, 
with a ~30% uncertainty. 

The value of the power-law index is similar to that obtained from C iv samples, a(C iv) 
~ 1.8, such as those drawn from the VLT-UVES Large Programme at z(Civ) = 2.16 
(Scannapicco et al. 2005) or from Keck- HI RES data, at higher redshift (Songaila 2001). 
The power-law fit of the latter (corrected for the different adopted cosmological param- 
eters) is also shown in figure 8. At N(Ovi) = N(Civ) = 10 13 5 cm~ 2 , the value of /(N) 
for O vi absorbers is nearly a factor of ten larger than that for C iv absorbers. 

4.2. Number density 

We use the power-law fit to /(N) for the O vi population to estimate the number density 
per unit z of Ovi absorbers: 




(4.1) 



dx = (i + z ) 2 {n A + n m (i + z) 3 }- - 5 dz, 



(4.2) 



or dX/dz = {(1 + z)/0.3} 5 when z > 1. 




(4.3) 
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Figure 8. Column density distribution of O VI absorbers. The dashed line is the fit to our data 
in the column density range 13.0 < log N(Ovi) < 14.3 (see text). The dotted line is the fit to 
the column density distribution of Civ absorbers given by Songaila (2001). 

We select conservative N(0 vi) limits of 10 13 and 10 15 cm" 2 , a range that is not drastically 
affected by incompleteness or sample variance. Adopting the fit with a = 1.71, we then 
get dn/dz = 74 at z — 2.3. Taking into account the range of possible values of the 
power-law index and normalization factor of /(N), we obtain 66 < dn/dz < 106. 

At low redshift, z = 0.1, surveys with FUSE and HST give dn/dz « 13 for a rest- 
equivalent width limit of w r ,min = 50 mA, or log N(0 Vl)=13.60 in the optically thin case 
(see Sembach et al. 2004). For this column density limit, we get dn/dz w 26 at ~z — 2.3, 
whereas we expect a somewhat higher value, dn/dz w 36, in the case of an uncvolving 
Ovi population. However, comparison between the values of dn/dz at z = 0.1 and 2.3 
is not straightforward as O vi absorbers may trace different populations at low and high 
redshift. 

4.3. Cosmic density o/Ovi absorbers 

The O VI cosmic density can be expressed as a mass fraction relative to the critical density, 
Pcrit- It can be estimated either from the individual, observed Ovi column densities or 
using the power-law fit to /(N) of Ovi absorbers to correct for incompleteness. 

4.3.1. Observed Ovi cosmic density 

The mean cosmic density of a given ion can be expressed as: 

ftb,ion - {H m ion /cp crtt }{J2 Nion/ AX ^ = 2 - 20 X 1(r22 £ N -n/ AX ^ ( 4 " 4 ) 

where Ho is the Hubble constant, nii on and Nion the atomic mass and the sum of the 
column densities of the given ion, respectively, and AA the total redshift path. For our 
O VI sample, we obtain f2t>(0 vi)= 1.51 x 10 -7 , a value higher than previous estimates by 
a factor 1.3 (Simcoe et al. 2004: sample restricted to O vi systems with strong, associated 
Hi absorption [see § 2.1]) and 1.8 (Carswell et al. 2002: two sightlincs, none with very 
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strong Ovi absorbers). The contribution of the Ovi type 1 population to fib(Ovi) is 
35%. 

4.3.2. Ovi cosmic density corrected for incompleteness 

The mean cosmic density of O vi ions can also be written as: 

O b = 2.20 x 10~ 22 J N/(N)dN (4.5) 

Using our fit with a — 1.71 and the same N(Ovi) limits as in § 4.2, we get fib(Ovi) 
w 3.5 x 10~ 7 , thus an incompleteness correction factor of 2.3 at z = 2.3. The uncertainty in 
the power-law fit of /(N) leads to values in the range 2.6 x 10~ 7 < ftb(0 vi) < 6.7 x 10~ 7 . 

To estimate the mean cosmic density of oxygen, fib(O), we use the O VI mean ionization 
level obtained in the pure photoionization case, (Ovi/O) = 0.15 (see § 3.2). Under the 
other ionization conditions investigated in § 3.3, this ratio is either similar or smaller. 
We then get fi b (0) « (2.3±§;J) x 10~ 6 . 

Using the solar oxygen abundance given by Anders & Grevesse (1989), yields: 

log (fi b (O)/fi b (O) ) = -2.22. 
This result demands attention for the two following reasons: (1) the above value is close 
to the median of [O/H] found for the Ovi type absorbers (IGM), [O/H] = -2.07, but 
well below that of the Ovi type 1 absorbers (metal-enriched sites), [O/H] = —0.33 (see 
figures 6 and 7), and (2) it is smaller than the mean metal enrichment of the IGM by 
star-forming galaxies at z ~ 2.5 (Z SF ) by a factor of about 3.7 and 6.6 when adopting the 
values of Z SF « 1/45 and 1/25 solar as given by Ferrara et al. (2005) and Pcttini (1999), 
respectively. Consequently, there is still a shortfall of observed metals as compared to 
those produced by LBGs, but about a factor three smaller than previously thought. 

Our sample contains very few cases of unambiguously broad Ovi doublets (b > 16 km 
s -1 ) which could trace hotter parts (T > 2.5 x 10 5 K) of metal-rich sites. However, if 
most of the gas in these sites is at even higher temperatures, oxygen will then mainly be 
in the form of O VII and O vm ions and not detectable with present-day X-ray satellites. 



5. Origin of the O vi absorbers 

5.1. Nearest strong Hi absorption system 

From a pixel analysis of VLT-UVES Large Programme quasar spectra, Aracil et al. 
(2004) found that weak Ovi absorption associated with weak Hi absorption (0.2 < 
r(Ly-a) < 1 or 12.9 < log N(Hl) < 13.6 for 6(Hl) = 30 km s" 1 ) is predominantly 
detected in the vicinity (Av < 300 km s _1 ) of strong Hi absorption (r(Ly-a) > 4). 
These authors suggested that the Ovi absorption arising in regions spatially close to 
strong Ly-a absorption may be part of outflows from overdense regions. 

The Ovi type 1 population from our study should exhibit the same property as the 
weak O vi absorptions analyzed with the pixel analysis method, since there is an overlap 
in their N(H i) range. The distribution of Aw between O vi or C iv-only type 1 systems 
and the nearest strong Ly-a system (r(Ly-a) > 4) is presented in figure 9. For the few 
cases of an O vi doublet associated with a saturated Ly-a line, Av is nul. It should be 
noted that all C iv-only type 1 systems have unsaturated, associated Ly-a lines. 

Among the Ovi and Civ-only type 1 systems, 64% and 63%, respectively, have a 
strong Ly-a system at Av < 450 km s _1 . Results from both the pixel analysis method 
of weak Ovi systems (log r(Ovi) ~ —1.35) and the study of individual Ovi absorbers 
suggest a link to gas outflows. 
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Figure 9. Distribution of the velocity difference between type 1 systems and the nearest 
strong Ly-a system, with the totals outlined. 



5.2. Gas density of the Ovi absorbers 

The gas overdensity of the O vi absorbers is estimated for two cases: photoionization by 
a hard UV metagalactic flux and hydrostatic equilibrium (Schaye 2001). 

In the photoionization case U is fixed by the O vi/C iv ionic ratio, assuming a relative 
O/C solar abundance. In the range 2.0 < z < 2.5, the adopted hydrogen photoionization 
rate is r(Hi) s=a 1.5 x 10~ 12 s _1 . Using the mean baryon density at each z(Ovi), we get: 

i(C/)!*4.0 [/-'([l + zl/3)- 3 . (5.1) 

The results are shown in figure 10 for the different types of O vi absorbers. The median 
values of 5(U) for the type (metal-poor) and type 1 (metal-rich) populations are equal, 
5(U) = 22, and that of the type 2 population is ^40% smaller. A Kolmogorov-Smirnov 
test shows that the types and 1 populations have the same 5(U) distribution at the 
97% confidence level. 

For the hydrostatic equilibrium case (Schaye 2001: equation (8)), we assume a gas 
temperature T = 4 x 10 4 K and the same photoionization rate as above. This gives: 

5(G) = 3.7 x 10~ 9 N(Hi) 2/3 ([l + z}/3)- 3 . (5.2) 

The results are presented in figure 11. Contrary to the photoionization case, there is 
a marked difference between the types and 1 populations. The median values of 5(G) 
are 41 and 4.6 for the type and type 1 absorbers, respectively, and that for the type 
2 absorbers is 8.3. Moreover, the values of 5(G) for ^80% of the type 1 population (log 
5(G) < 1.1) do not overlap with those obtained for the type population. 

We now compare in figure 12 the values of 5(U) and 5(G) to check the validity of 
the assumption of hydrostatic equilibrium. For the type absorbers, a Spearman rank 
correlation test shows that 5(G) and 5(U) are correlated at a >99% confidence level. The 
mean value of their 5(G) /5(U) overdensity ratio is close to 2.0. It cannot be substantially 
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Figure 10. Distribution of the O VI absorber overdensity in the photoionization case. 
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Figure 11. Distribution of the Ovi absorber overdensity in the hydrostatic equilibrium case. 



decreased as its dependence on T(Hi), mass fraction in gas and T is small (8(G) / 6(U) oc 
T^TCRi)- 1 / 3 (O b /^m)" 1/3 )- Therefore, this departure of 5(G)/5(U) from unity may 
suggest that a fraction of the observed H I is not in the O vi phase. Nevertheless, the 
correlation between 8(G) and 6(U) suggests that the O vi type population is roughly in 
hydrostatic equilibrium. For the type 1 absorbers, 6(G) and 6(U) are totally uncorrelated 
which implies that hydrostatic equilibrium is not a valid assumption for this population: 
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low H I column density absorbers do not trace low density regions of the IGM. This further 
supports that the metal-rich and metal-poor absorbers trace different populations. 




1.0 1.5 
log 5(U) 



Figure 12. Comparison of the O VI absorber overdensities in the photoionization versus 

hydrostatic equilibrium cases. 



6. Conclusions and prospectives 

Our large VLT-UVES sample of 136 O VI absorbers at z — 2.28 towards ten quasars 
enables a study of the highly ionized phase of the IGM, in particular its metal enrichment 
and contribution to the cosmic baryon density. 

Previous Ovi studies at high z have uncovered a few systems with high [O/H] (> 
— 1.0) abundances, motivating us to fully investigate this class of metal- rich absorbers. 
Because those systems already identified have low Hi column densities (Bergeron et al. 
2002; Carswell et al. 2002), our sample includes all detected Ovi systems whatever the 
strength of their associated H I absorption. In contrast, the survey by Simcoe et al. (2004) 
only includes systems with N(Hi) > 10 13 6 cxnT 2 . We restrict our sample to absorbers 
with both lines of the Ovi doublet clearly detected or, if partially blended, enough 
unambiguous structure to allow for dcblcnding from Lyman lines. 

Since nearly half of the Ovi absorbers have small line widths, b < 12 km s _1 or 
T < 1.4 x 10 5 K, photoionization must be the dominant ionization process. We thus 
introduce an observational identification critcrium to separate the classes of metal-poor 
(type 0) and metal-rich (type 1) absorbers. Selecting a hard UV background flux (see 
Bergeron et al. 2002; Carswell et al. 2002 ; Tclfcr et al. 2002) and assuming a 0.1 solar 
metallicity yields a column density ratio N(Ovi)/N(Civ) < 0.25 and > 0.25 for the type 
and type 1 populations, respectively. 
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The bulk of the b distributions of these two Ovi populations are similar except at 
the highest velocities. However, we stress that very few individual components are un- 
ambiguously broad, a result of blending for complex, multiple systems and limited S/N 
(~ 30-40) in the O vi range. 

The cosmic oxygen abundance is derived under the assumptions of photoionization, 
coexistence of O VI and C IV in the same phase, and a solar O/C relative abundance. The 
overall [O/H] distribution is clearly bimodal with median values of [O/H] equal to —2.05 
and —0.33 for the type and type 1 populations, respectively. This is not a consequence 
of a strong difference in ionization levels between the two types of O vi absorbers. All 
of the type 1 O vi systems and all but two of the type systems have associated C iv 
absorption (detection limit of N(Civ) « 1 x 10 12 cur 2 ). Their N(0 Vl)/N(C IV) distri- 
butions, which cover about two orders of magnitude, are similar with median values both 
close to 10. A high metallicity (median [O/H] > —1.0) for the type 1 population is still 
found under different ionization conditions: photoionization together with either a gas 
phase of constant density (overdensity 6 = 10 at z = 2.2) or a temperature fixed by the 
line width of the main Ovi component of each system (always < 14 km s _1 ). 

The N(N v) /N(0 vi) ratio cannot be used to constrain the ionization level of the O vi 
phase because the N/O relative abundance departs from the solar value. In most Ovi 
absorbers, the Nv doublet is weak or undetected and the nitrogen abundance relative to 
oxygen is usually well below solar (Bergeron et al. 2002). For a very few Ovi systems, 
the strength of the N v absorption is similar to those of C iv and O vi and the nitrogen 
abundance may be enhanced relative to that of oxygen (Carswell et al. 2002), as also 
observed in quasar associated systems (e.g. Hamann et al. 2000). In our type 1 Ovi 
sample, associated Nv absorption is either absent or very weak, except in one case 
already reported by Bergeron et al. (2002). This absorber at z — 2.352 in Q 0329—385 
was labelled as "intrinsic" by these authors because its properties are typical of those of 
associated systems, even though it is at 6200 km s _1 from the quasar emission redshift. 

Our O vi sample is large enough to derive the first estimate of the O vi column density 
distribution, although incompleteness becomes evident at N(Ovi) < 1 x 10 13 cm~ 2 and 
sample variance may be important at N(0 vi) >2x 10 14 cm~ 2 . A power-law fit, /(N) oc 
N~ Q , yields a s=a 1.7 ± 0.5, a value similar to that found for Civ samples, a(Civ) ~ 1.8 
(Songaila 2001; Scannapieco et al. 2005). In contrast, the normalization factor, /(N) 
= 2.3 x 10" 13 (with an uncertainty of - 30%) at N(Ovi) = 10 13 5 cm~ 2 , is about ten 
times larger than that of C iv absorbers. We aim to better constrain /(N), particularly at 
large O vi column densities, by analyzing a larger number of sightlines in a future paper. 
There we will also include blended O vi components associated with C IV absorption. 

We use the fit to /(N) in a conservative N(Ovi) range, 10 13 -10 15 cm~ 2 , to estimate 
the number density per unit z of O vi absorbers as well as their cosmic density. This is 
a first step for correcting dn/dz and fib(Ovi) from incompleteness and sample variance. 
We find dn/dz = 74±| 2 . Selecting an integration range as wide as is usually adopted 
for Civ (10 12 -10 16 cm~ 2 ) would yield a larger value of dn/dz (and of fi b (Ovi)) but 
the uncertainty on the result would then be far too large. In the case of an unevolving 
population and a column density lower limit equal to that of low z = 0.1 Ovi surveys, 
that is N(Ovi) min = 10 13 6 cm~ 2 or w r , m in = 50 mA (see e.g. Sembach et al. 2004), the 
expected value of dn/dz at z — 2.3 derived from the low z samples is 36 (assuming a 0.1 
solar metallicity and an ionic fraction Ovi/O = 0.2), whereas that obtained from our 
/(N) distribution equals 26. The suggested decline of dn/dz with z is not straightforward 
to interpret as the Ovi absorbers may trace different populations at low and high z. 

The Ovi cosmic density estimated from the individual, observed column densities is 
Ob(Ovi)Ri 1.5 x 10~ 7 , i.e. higher than previous estimates by a factor 1.3 (Simcoe et al. 
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2004) and 1.8 (Carswell et al. 2002). This increase is due to the high contribution (35%) 
of the type 1 population to f2b(0 Vl). The value derived from the /(N) distribution is 2.3 
times larger: f2b(Ovi)« (3.5±q |) x 10~ 7 . This illustrates the effects of incompleteness 
and sample variance in our sample, even within the conservative N(0 vi) range adopted. 
To get the element cosmic density, we use the mean ionic fraction obtained in the pure 
photoionization case, (Ovi/O) = 0.15, which yields ^b(O) ss (2.3±g g) x 10~ 6 . 

Adopting the solar oxygen abundance given by Anders & Grevesse (1989), we get 
log (£!b(O)/£!b(O)0) = —2.22. This value is well below that of the metal-rich population 
and also smaller than the metal enrichment of the IGM expected from high z star-forming 

galaxies, ([O/H]) 1.40 (Pettini 1999) or -1.65 (Ferrara et al. 2005). Although the 

problem of missing metals at high z (where previously an order of magnitude disparity 
was measured) is now less severe as a result of our O vi survey, there remains a shortfall of 
observed metals by about a factor of four as compared to those produced by star-forming 
galaxies. 

Other properties of the type 1 Ovi absorbers suggest a tight link to galactic halos. 
This population is predominantly detected in the vicinity (Av < 450 km s _1 ) of strong 
Hi systems (r(Ly-a) > 4). This is also the case for Civ-only metal-rich absorbers (Ovi 
doublet cither outside the observing range, z < 2.0, or fully blended with saturated 
Lyman lines). In the photoionization case, the type and type 1 Ovi absorbers have 
the same gas overdensity distribution, with a median value 5(U) = 22, but under the 
assumption of hydrostatic equilibrium the gas overdensity, 5(G), distributions of these 
two populations barely overlap. Moreover, the values of 5(U) and 5(G) are totally uncor- 
rected for the metal-rich population, whereas they are well correlated for the metal-poor 
population. Consequently, the assumption of hydrostatic equilibrium is not valid for the 
metal-rich O vi population: these absorbers do not trace low density regions of the IGM 
but rather gas outflows in the vicinity of active star-formation sites. 

If most of the gas in the metal-rich sites is at high temperature (T > 5 x 10 5 K), as 
suggested by Pettini (1999), oxygen will mainly be in the form of Ovu and Ovin ions 
and their signatures in the very soft X-ray range are not detectable with present-day 
X-ray satellites. For a phase at lower temperatures, 2 x 10 5 < T < 5 x 10 5 K, the Ovi 
and Hi species, but not Civ (ionic fraction Civ/C < 10~ 2 ), should be detectable. We 
have begun to search for these absorbers with broad (b > 50 km s _1 ), weak Ly-a lines 
associated with semi-broad, weak Ovi doublets (b > 15 km s _1 ). This is coupled to a 
statistical analysis of the Ly-a forest in simulated spectra (in progress). 

We also plan to acquire deep, multi-band images of the quasar fields with several 
metal-rich Ovi absorbers. If this population does indeed trace hot galactic halos, we 
expect to find a strong correlation with star-forming galaxies. Using these images together 
with spectroscopic follow-up of the associated galaxies may help clarify the ejection 
mcchanism(s) responsible for the metal-pollution of galactic halos and the surrounding 
IGM. 
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